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1.1 Overview and basic discussion of equations of stellar 

structure 

This section is intended merely as a reminder of things you have already 
learned (or will soon learn) in a first course in stellar structure and evolution. 
Here, we will introduce, and in some cases derive, the four basic equations of 
stellar structure. If you'd like additional details, by all means consult one of 
the standard texts in this subject - 1 recommend without hesitation the text by 



Hansen et al.l(l2004l) . 



1.1.1 Dependent and Independent Variables 

In solving the equations that describe the conditions within stars, we need to 
decide on a coordinate system. A reasonable simplification of the equations 
arises if we assume that stars are spherically symmetric. With that assumption, 
a single positional variable will suffice as the independent variable. Our choice 
could be r, the distance from the center of the star, or nir, the mass contained 
within a shell of radius r. For some purposes we may use r and for others, 
rrir- Their values are related through the equation of mass conservation (or 
continuity): 

^=47rrV(r), (1.1) 
dr 

where p{r) is the mass density at position r. 

As for the quantities that we would like to use to describe the physical con- 
ditions within the star at each position, we will define them as follows: 

• velocity: v (=0 for a hydrostatic star) 
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• density: p or n = Nap/ M 

• pressure: P 

• temperature: T 

• chemical composition (mass fraction): Xi 

• ion / charge balance: yi, rie 

• internal energy: U 

• entropy: S 

• heat flow parameters / cross sections: Krad, ^cond 

• energy flow: Li , Fconv 

• energy generation / loss: enuc, f^v 

These quantities provide a fairly complete description of the state of the 
stellar interior, but they clearly are not independent of one another. In fact, one 
can show that only four dependent variables at each point in a model suffice 
to fully describe the conditions there. For computation of stellar structure and 
evolution, those quantities are usually chosen as r (or m^), P, T, and Lr. 

Of course, the other quantities are important, but they can be derived or 
computed given the main four quantities. That said, the computation of a stellar 
model requires knowledge of the values of all of the above quantities. They 
appear in the coefficients of the equations, and contain all of the physics of 
relevance to stellar interiors. Thus they provide the direct connection between 
fundamental physics and the appearance and behavior of the stars. 

Their appearance in the coefficients of the equations of stellar structure al- 
lows us to further categorize what we need in terms of the primary mechanical 
quantities r (or rrir) and P, and the primary thermal quantities T and Lr. Nec- 
essary extra information is given by the composition {Xi). Various branches 
of physics provide the tools do determine important derived quantities include 
equation-of-state parameters p, /x, U, and S, from which we can calculate 
Vad> specific heats, and other thermodynamic quantities. From atomic physics 
we can compute yi and rig, and the opacities Krad and Kcond- Nuclear physics 
provides the tools for computing Cnuc and eneut> and hydrodynamics gives us 
hope for computing .Fconv 



1.1.2 The Equations of Stellar Structure 

Let's now see how the above set of parameters appears when we write down 
the equations of stellar structure. This is possible to do in a compact form only 
if we expect that the assumption of spherical symmetry everywhere is a good 
one (so that all quantities are functions only of r), and that the star that we 
are modeling is in hydrostatic equiUbrium (ensuring that time derivatives are 
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zero). Then and only then do the equations take on the form that is used, almost 
universally, used to produce stellar models for asteroseismological (and indeed 
most other) purposes. 



Mass Conservation 

The first equation has already been introduced above; equation 11.11 ensures 
that the mass and radius are consistent with the density from the center to 
the surface. Equation 1 1.1 1 also serves as a coordinate transformation between 
using rrir and r as independent variables. Consider a general quantity Z (which 
could stand for temperature, pressure, or something else). Then the differential 
equation describing the dependence of Z on r can be written down in terms of 
the dependence of Z on m,. as follows: 

dZ , y , . dZ 

— = Anr^p{r)— . (1.2) 
dr drrir 



Mechanical Equilibrium 

We're operating under the assumption of nothing moving... that is, there is me- 
chanical (hydrostatic) equilibrium within the entire stellar model. This means 
that the downward force of gravity at any position r, 

p(r)g(r) = p(r)^^ (1.3) 
must be balanced by the balance of the pressure upwards: 

P{r) - P{r + dr) = - (^^^ dr. (1.4) 

If they do indeed balance, then this equilibrium requires that 

dP Grrir , dP Grrir (15) 

dr r^ dnir 47rr^ 

Note, in passing, that equations ! 1 . 1 l and [T3] contain only r, m^, P, and p. So, 
under circumstances where one can write the pressure in terms of the density 
only (i.e. without an explicit temperature dependence) then these two equa- 
tions, plus the P{p) relation, suffice to completely describe the mechanical 
structure of the model. As an example, consider a polytropic equation of state 
like 

P(r) = Kp''{r). (1.6) 



Under this condition, then the two equations can be written in terms of the 
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independent variable (r) and two dependent variables (P and rrir): 
^^Grur 1 „ 

dr r2 if 1/7^ ^'^^ ^'-'^ 
^=W^PV.(,) (1.8) 

That is, two equations, two unknowns - so just add a pair of boundary con- 
ditions, and you have a complete stellar model, in hydrostatic equilibrium! 

Energy Generation 

In the more common case where we do care about the thermal content of stellar 
material, we need to determine the energy balance within each zone of the star 
To this end, we appeal to conservation of energy: that is, the energy flowing 
into a zone must be balanced by the energy flowing out of that same zone, 
possibly affected by energy production or consumption within the zone. If no 
energy is gained or lost, then clearly 

(1.9) 

But more generally, 

Lr+dr = +47rrVW ("^ ^r, (1.10) 

where e is the net energy generation rate per unit mass, and Q is the heat con- 
tent per unit mass. The quantity e generally denotes nuclear processes (i.e. 
energy production via fusion) but also includes energy losses through neutrino 
emission. In those cases, e can be computed as a function of the thermody- 
namic state and the composition, i.e. e(p, T, X^). 

The term involving the time derivative of Q describes the rate of heat gain 
or loss, per gram, of material. Generally, we have 

dQ^dE^ dV^dE Pdp 

dt dt dt dt dt ' ' 

where the first term on the right is the rate of change of internal energy, and the 
second term accounts for any PdV work done on (or by) the zone. Equation 
ll.lll can be recast using the definition of entropy, 

which leads to the more familiar form for the energy conservation equation 
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We note that the time derivative in equation ll.l3l is the only place where time 
explicitly appears in the equations of stellar structure. More on that later 



Energy Transport 

Finally, we require an equation that tells us how temperature changes with 
position in a stellar model. To get there, we invoke thermal equilibrium - the 
assumption that since energy is flowing out of the "top" of a given region, 
energy must also be flowing in from the bottom. We model this transport as 
a diffusive process driven by the fact that the energy density has to change 
through the region. This in turn is a result of the continuously increasing radius 
of each zone along with any temperature change. Schematically, we can write 
the radiant flux Fr as follows: 

Fr = energy-density gradient x speed x mean-free-path. 

In a less schematic form, we have 



dr 



xcxA (1.14) 



where A is the mean free path between photon scatterings, and c is the speed 
of light. Recognizing that for radiation, the mean free path A = (np)^^ where 
K is the opacity per gram, we can rewrite the above as 

4ac q dT 

Fr^—T'— (1.15) 
Snp dr 

for the energy flux carried by radiation. Multiplying by the surface area of the 
zone and rearranging yields an expression for the gradient in the temperature 
when radiative diffusion carries the flux: 

dT ^ 3kp Lr 
dr 16nacr^ 

We can also express the temperature gradient in more general terms as a 
function of the pressure gradient: 

dT dT dP T d\-aT dP 



If we define 



then 



dr dP dr PdlnPdr 
_ dlnT 



d\nP 



(1.17) 



(1.18) 



^^-V^f. (1.19) 

dr r^ P 



6 



This is a general expression for the temperature gradient because the mecha- 
nism of heat transport is not specified, but is contained in the way V is calcu- 
lated. 

Equation 1 1.161 in fact, can easily be transformed to look like equation ll.191 
Looking at equations 11.161 and 11.191 we can define a "del" for the case when 
radiation carries the flux: 

^rad = -1^ '^TTJT- (1-20) 

loTrac 1 * uMr 

In turn, we can now write the temperature gradient in the model as 

dT _ GMrpT 

~r — ^Vrad 5 fT (1-^1) 

for the case where V is determined by radiative diffusion (that is, V = Vrad)- 
Later, we will consider other forms of heat transport which can provide values 
for V under a variety of physical conditions (for example, when the material 
is convective or conductive). 



1.1.3 The Constitutive Relations - Where the Physics Is 

The previous section summarized the four differential equations of stellar struc- 
ture and the general background from where the came. Each equation relates 
the dependent variables to the independent variable, but each also includes 
some other factors and terms, within which the physics that governs stellar 
structure are reside. In this section, we outline these quantities and describe 
the "constitutive relations" that provide the route to evaluating the values of 
these terms. 

The Equation of State 

The dependent variables P and T, along with the compositional mix of the stel- 
lar material (mass fractions Xi for each atomic species i) suffice to determine 
the density p, the ionization state of the material the relevant thermal quanti- 
ties such as internal energy, and the ionization state of the material. The route 
to calculating these quantities can be a difficult one in complete generality, 
but with a few simplifying assumptions we can make progress. The principal 
assumption is that the material in a stellar interior is everywhere in local ther- 
modynamic equilibrium, and that all quantities are isotropic. The fact that the 
photon mean free path is much smaller than the length scales of interest (the 
pressure scale height, stellar radius, etc.) ensures that this approximation is a 
relevant one. 
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With that assumption (and assuming isotropy) one can show that for an 
unionized gas, the pressure and density are related through 

P^Yl "^'^^ " NapX,/A^. (1.22) 

i 

We now define the "mean molecular weight" /i as 

/i"' = ^X»M, (1.23) 

i 

SO that equation ll.22l becomes 

P^P-NAkT (1.24) 

and also the internal energy per gram is 

E^\-. (1.25) 
2 p 

Of course, the interiors of stars are mostly ionized, so that the above equa- 
tions need some modification to account for the free electrons in addition to 
the nuclei. Thus we have 

Pe + Pl = Pe + —NAkT (1.26) 
A*/ 

and further note that the right-hand term for Pj remains valid even in the case 
of electron degeneracy when dealing with most stars that have asteroseismic 
potential. On the other hand, for the electron pressure under non-degenerate 
conditions, we have 

P, = —NAkT (1.27) 

where pe depends on the ionization state of the material. 

As the reader might expect, there are complications beyond simply comput- 
ing the ionization state. Departures from ideal gas behavior can be significant 
under conditions found within the Sun and the stars, and we will discuss those 
later. In massive stars, radiation pressure (Prad = becomes increas- 

ingly important in the outer layers as the mass increases. 

Energy Generation Rates 

Stars are, through much of their visible Ufetime, gigantic nuclear furnaces. To 
ensure energy conservation as per equation 11.131 we need an accurate deter- 
mination of the net energy generation rate per gram, e. For the nuclear fusion 
process, this means computing the rate of interaction between a target and a 
projectile. We express the probability of these interactions in terms of a cross 



8 



section a for a given process. Then the rate of the interaction is proportional to 
nav, where n is the number density of targets and v is the relative velocity of 
the reactants. The velocity will depend on energy (and scale with \/fcT). 

The reaction cross section cr depends on energy as well, and theory can pro- 
vide a reasonable functional form of that energy dependence, but the scaling of 
that relation requires laboratory measurement. Since the energy dependence of 
the velocity and cross section can be understood, conmion practice is to tabu- 
late ((Tw), averaged over the energy distribution corresponding to a temperature 
T. 

For resonant reactions, your stellar structure class should have taught you 
that 

{av),^ = K, g^T-^/^e-^^/T (1.28) 

where Ki and K2 are constants that depend on the properties of the interacting 
particles, and the lifetimes of the incoming state and outgoing state are indi- 
cated by the energy widths Fj and Tj, with the total hfetime given by l/F. The 
quantity (/ is a statistical weight factor. In the case where a reaction does not 
have a resonance in the energy range of interest, the {av) value is given instead 
by 

Here, 5(0) is the "astrophysical S factor" which is the value, at zero energy, 
of a slowly varying function of energy that helps isolate a part of the energy 
dependence of the reaction. Once we have determined these cross sections, the 
rate of energy production per gram of material can be computed with knowl- 
edge of the energy yield per reaction and the abundances of projectiles and 
targets: 

ajip, T, X) = Qij pNl^ ^{^v)i^ (1.30) 

Radiation, Conduction, and Convection 
Radiative Transport 

Section 1.2.4 expressed the temperature gradient in terms of V, which is 
the local power law slope of the T{P) relation. The value of V will reflect 
the dominant energy transport mechanism. The example shown in that sec- 
tion was for energy transport by photon diffusion, which defined a Vrad- That 
quantity depends in turn on the conditions at that point (P, T, M,., L), some 
constants, and the radiative opacity Kr{p, T, Xi). The radiative opacity repre- 
sents scattering photons through a variety of atomic processes involving elec- 
tronic transitions (bound-free, bound-bond) and free particles (free-free and 
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electron scattering). Though the scattering cross-section will be a sharp func- 
tion of frequency (we live in a quantum world) the radiative opacity to be used 
in computations of isotropic stellar interiors may be simplified using a suitable 
average opacity over the expected frequency or wavelength distribution of the 
flux at the local temperature T. A favorite frequency average is the Rosseland 
mean opacity 

-L _ Jo K„ dT "'^ 



Jo dT "'^ 



(1.31) 



where the average is weighted by the local flux derivative (and suitably nor- 
malized). Note that this is an average of the reciprocal of the opacity, since the 
radiative flux is inversely proportional to the opacity (see equation 1 1.151 1 and it 
is F that is the quantity that is driving the temperature gradient. 

Conductive Transport 

Under certain circumstances in stellar matter (as well as the much more com- 
mon household applications), heat can be transported by conduction. Under 
the constraints of the formalism adopted in Section 1.2.4, we must seek a con- 
ductive opacity Kc that encapsulates the physics of conduction. Thus, we write 

3pDe 

where the diffusion coefficient for conductive heat transport is 

D,^^. (1.33) 

In equation 11.331 Cy is the specific heat at constant density, Ve is the thermal 
velocity of the electrons, and A is again a mean free path. Skipping several 
steps, we can eventually show that 



cx !^Zf - (1.34) 
M/ \Pj 

where the constant of proportionahty is the hard part. Finally, we note that 
when both conductive and radiative opacities are comparable, the combined 
opacity adds in parallel; that is since 

,.,,^r^.-^J^iL(L^L\ ,,35) 



3p dr \ Kr 



the net opacity is 



(1.36) 
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Convection 

Finally, we come to one of the continuing perplexities of computational stel- 
lar evolution, and that is the treatment of the convective flux. Again, as the 
energy transport problem has been posed, the problem reduces to determining 
the value of V when convection carries the flux - that is, when heat transport 
via bulk turbulent motion, we employ Vconv in equation |l.l9l 

Convection is required should the radiative temperature gradient exceed the 
adiabatic temperature gradient. To demonstrate this, consider an ephemeral 
blob of material where the temperature exceeds the ambient surrounding tem- 
perature. Such a blob will have a lower density than the surroundings since 
pressure equilibrium is enforced. Archimedes taught us that material with a 
lower density than its surroundings experiences a buoyancy force upwards; 
this blob begins to rise. If it does so without losing its identity by exchanging 
heat with its surroundings, then the temperature will drop as it moves upwards 
to lower pressures. The temperature will fall along an adiabat. The surround- 
ings, however, need not be along the same adiabat; the degree to which the 
surrounding temperature falls is determined by the local value of V and is pre- 
sumably Vr. If, after a displacement upwards, the blob finds that it is warmer 
than its surroundings, it will continue to experience an upward force. That 
can only happen if the surrounding temperature drops faster than the adiabatic 
(blob) value - this occurs when V ad < V^. If not, and the blob is cooler than 
its surrounding after an upward displacement (Vad > V^), then the blob's 
density will become higher than the surroundings, and it will sink back to its 
equilibrium value. 

Eventually, a blob that is buoyantly unstable (and therefore has a higher 
heat content than its surroundings) will eventually dissolve, resulting in a net 
transport of heat from the inside of the star to the outside. This excess heat 
transport will lower the overall gradient somewhat; a cascade of such blobs 
therefore can globally affect the temperature gradient in the star . . if Vr is 
larger than V^d to start. 

This convective instability, therefore, is generated in regions of the star 
where is large when compared to the adiabatic value - from eauation ll.201 
that is, where k is large or where L/rrir is large. In the first case (large Kr) we 
often see the effects of convection in cooler regions of stellar envelopes where 
partial ionization of a dominant species can drive up the radiative opacity. In 
the latter case, if energy generation is confined to a small region, so that the lu- 
minosity rises rapidly over a small mass range, then convection can be required 
to carry all of this locally generated flux. This is the situation in the cores of 
stars, and in helium shell burning on the asymptotic giant branch. 

In the extreme, convection suppresses the temperature gradient in unstable 
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regions, bringing it down to Vq^ in the case of perfectly efficient convection. 
Of course, modeling this process as resulting in neutrally stratified material is 
a simplification of the turbulent and dynamic motions expected where convec- 
tion is present. Approximating Vconv as Vad in regions that are unstable, and 
mixing the material completely to simulate these dynamical processes, may be 
good enough for some purposes (such as core convection). And we have the 
added benefit that Vad is a thermodynamic quantity, meaning that it comes for 
"free" from the equation of state constitutive relation. 

But in reality we should allow for imperfect convection, and seek a truer 
treatment of convection that can yield a value for V that is intermediate be- 
tween Vr and Wad- One such treatment, with a small number of free parame- 
ters, is the famous mixing-length theory. Tuning the mixing length to pressure 
scale height ratio (a = £/Xp) so that models reproduce the solar radius and 
luminosity at the solar age, for example, provides a touchstone to explore how 
convection might work in other stars. Many variants of the mixing length have 
been proposed and employed over th e 52-H years since it w as introduced into 
astrophysics by Erika Bohm-Vitense (IBohm- Vitensel 1 1 95 8h . 



1.1.4 Conclusions - room for improvement? 

Given the importance of convective energy transport across stellar astrophysics, 
and the relatively crude model we have been using for it for decades, it is not 
surprising that astrophysics is being held back by our ignorance of the convec- 
tion process. Computational hydrodynamics simulations of conditions found 
in stars are being brought to bear on the subject, but are (still) computational 
expensive and time consuming. Implementation into stellar evolution codes is 
a relatively young endeavor, but an important one that is relying on seismic 
diagnostics of the Sun and stars. 

As we will see in the remaining portions of my contribution, the coefficients 
of the equations of stellar structure rely on basic physics that, in some cases, is 
well understood. However, like convection, many of the processes that go into 
these coefficients are not very tightly constrained by experiment or by theory. 
With asteroseismic probes, we may be able to place new and strong constraints 
on a variety of physics problems. 



1.2 The physics behind the coefficients 



In the previous section, we reviewed how the basic equations of stellar struc- 
ture result from various aspects of equilibrium in a self gravitating (mostly) 
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gaseous sphere. The coefficients within those equations connect with some ba- 
sic physics: the equation of state of material under extreme conditions, nuclear 
reaction cross sections, and the interaction between matter and radiation. In 
this section we review the relevant constitutive relations and point out areas 
where the physical inputs could be improved through the indirect probes pro- 
vided by asteroseismology. 



1.2.1 The Equation of State: T, Xi) 
Basic elements 

The "easy" form of the equation of state involves the assumption of a perfect 
gas (perhaps undergoing ionization). Complicating factors include the role of 
electron degeneracy, but the computational aspects of that are well known. Less 
manageable aspects include non-ideal effects, where some (electromagnetic) 
interaction between particles leads to effects like pressure ionization. Coulomb 
interactions, and state changes such as crystallization. 

To begin with, though, the ideal gas relation leads to an equation of state of 
the form shown in equations 11.261 and 11.271 where He depends on the ioniza- 
tion state of the material. Ionization balance, in the ideal case, is provided by 
recursive solution of the Saha equation 

nk+i^^ = ^k+i{T) 2(27rfcTme)3/2 ^_^^^(fer^^ ^^^^^^ 
rik " Uk{T) 

where the index k denotes the ionization level, Xk is the ionization potential of 
level k, and Uk is a suitable partition function. 

These all work well enough when the density is "sufficiently" low. But as 
the density increases, the quantum statistical nature of the electrons begins to 
overwhelm the simple ideal, perfect gas assumption. For Fermi-Dirac particles 
such as electrons (with available spins ±1/2), the number density is a function 
of particle momentum: 

2 1 

h-^ exp + mc^ + h[p)\/ki + 1 

where, generally, 

E{p) ^ {p^c^ + m^c^ f^ - mc^ (1.39) 

and the chemical potential is /i. 

From these, one can obtain the pressure and internal energy from 

P=i 1 n{p)pv4:TTp'^dp -a.nd-E^ j n{p)E{p)4:TTp'^dp . (1.40) 
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The interesting part of the above is the energy (and temperature) dependence 
(the left-had term of equation ll.381 which allows us to write 

n = — / F{E)p^dp (1.41) 
Jo 

If we now define the Fermi energy in terms of the Fermi potential '■ 

Ep = -mc^ (1.42) 

then in the limiting case of low T (that is, kT < Ef) is F{E) = l\i E < 
Ep and F{E) — \f E > Ep- At these low temperatures, then, P and E 
become independent of T and depend only on the Fermi energy (which in turn 
depends only on the density). This is the well known degeneracy pressure for 
a r = system. At higher temperatures, thermal energy can bump particles to 
E > Ep at the top of the energy range only, eroding the "edge" of the F{E) 
distribution at the E = Ep edge. Computation of the electron pressure under 
this intermediate-degeneracy range is well understood (but does involve some 
messier integrals). 



IVouble in an ideal world - charged particles do interact 

In equation 11.371 the relative populations depend on the ionization energies 
Xfc. However, those energies are valid for the ion in isolation, and are mea- 
sured relative to the continuum. At densities relevant to stellar interiors, how- 
ever, the electric field from surrounding ions can influence the ionization level 
by effectively suppressing the continuum, making ionization easier at a given 
temperature. The separation between ions, a, depends on the density: 

a = 7.3 X lO^^cm X f p— j (1.43) 
so that the surrounding ions lead to a depression of the continuum of 

Ax« ^ =9.8eVx Z2 U-j . (1.44) 

For hydrogen, at densities of a few g/cm^, the effective change in xh is com- 
parable to the ionization energy itself - that is, hydrogen should be completely 
ionized independent of the temperature. It is this "pressure ionization" that is 
responsible for the interiors of stars being fully ionized. Without it, the high 
densities within stellar interiors could lead to recombination if the Saha equa- 
tion remained unmodified - that is, if the unadulterated ionization potential is 
used. Another consequence of these proximity effects is that the internal en- 
ergy is reduced below perfect gas levels. Simply, the ion-electron interactions 
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result in a binding energy that would not exist otherwise, adding a negative 
component to the total internal energy. 

If the ions and electrons "feel" one another, then at some level, there are ion- 
ion interactions as well. This is because the shielding of ions by surrounding 
electrons can become imperfect if the density is high enough. Under those 
circumstances, we must consider Coulomb repulsion between the imperfectly 
shielded ions. The energy of these Coulomb effects be proportional to Z^e^/a. 
Scaling this to kT gives us the Coulomb parameter Tc- 

Vc = —-=2.21Z^4k] TTT^ TTT (1-45) 



akT ■ Vl06/ VlO^y \12, 

using cgs units for scaling the density, Kelvin for the temperature, and AMU 
for the atomic weight A. When Fc = 1, Coulomb effects begin to be signif- 
icant; in the solar interior, Tc = 0.1. While still small, the sensitivity of he- 
lioseismic analysis demands that Coulomb effects be taken into account. For 
other stars, such as red giant cores and white dwarfs, Tc can be much greater 
than 1 , and mutual ion repulsion becomes an important element of the equation 
of state budget. 

In particular, as first proposed bv ISalpetei (1961), in white dwarfs where 



Tc can exceed 100, this mutual ion repulsion could result in crystallization of 
white dwarf interiors. In the modern formulation, the onset of crystallization 
is believed to occur when Tc exceeds 175. In that case, when a white dwarf 
interior cools to below a threshold temperature 

T^tai « 3.4 X 10^ V f ( TTTF^") K (1.46) 



16 J V 10*^8 cm 

the core solidifies. Note that the crystallization varies as with Z^/A^^^^, mean- 
ing that as a white dwarf cools, heavier elements crystallize at higher temper- 
atures, so in the course of white dwarf cooling, they crystallize before lighter 
elements. The order of crystallization remains a difficult issue, because of ro- 
tation and other mass motions in the degenerate interior. Whether the core of 
the white dwarf solidifies as some sort of alloy (principally of oxygen and car- 
bon) or if the oxygen precipitates out prior to reaching the lower temperature 
of carbon crystallization, remains an interesting and difficult question. 



Full-up Equation of State Calculations 

Figure 11.11 summarizes the various places in phase space where these effects 
become important. Putting together the elements of Sections 11.2.11 and 11.2.11 
the computation of a complete, realistic equation of state is not trivial. One 
must address ionization balance using a full mixture of elements and realistic 
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Figure 1. 1 Equation of state regions in the p — T plane, adapted from 
iPaxton et al. (2011). The labeled lines show borders where various effects 
are relevant as described in the text. The sold lines also denote the run of 
density and temperature for stellar models on the main sequence and for a 
O.8M0 white dwarf. 



partition functions for the bound states, while also accounting for degener- 
acy and non-ideal effects. Most modern treatments use variations of the free- 
energy minimization approach to compute extensive tables with equation of 
state parameters for mixtures of interest at a range of pressures (or densities) 
and temperatures. 

As an example, the OPAL effort (see ( Rogers et al. , 19961) and references 
therein) covers conditions relevant to the Sun and solar -type stars. For lower 
temperatures and relevant densities. lSaumon et al.l d 1995b produced widely used 
tables. Under conditions relevant to white dwarfs and red giant cores, the 
most recent useful algorithm (with analytic fits to the full calculations) is by 



Potekhin and Chabrieii (1201 Oh which allows for Coulomb interactions as well 



as crystallization. 
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Figure 1.2 Astrophysical S factor for a representative nu- 
clear reaction as measured. Data from the NACRE compilation 
(http://pntpm3.ulb.ac.be/Nacre/barre_database.htm) for the ^^C(p, 7)^'^A^ 
reaction. Typical relevant stellar energies are approximately 1-2 keV, well 
below the experimentally accessible energies. 



1.2.2 Energy Generation and Loss Rates: e{p, T, Xi) 

Nuclear cross sections 

Developing expressions for the rate of energy generation through nuclear fu- 
sion (and loss through neutrino emission) is at the crossroads between exper- 
imental and theoretical nuclear physics and computational astrophysics. The 
essential scheme is a collision-physics one: this is a "rate = n{(jv)" problem 
at its core. The relative velocity v depends on the mean kinetic energy per 
particle {v{E) = yf{kTlm)) assuming a Boltzmann distribution for particle 
velocities, and a{E) is the "cross section" for the given interaction. 

Nuclear physics provides a functional form for the cross section a{E), but 
laboratory measurements are essential for an absolute determination for use in 
stellar models. Generally, these semi-empirical cross sections are tabulated as 
{<7v)ij for interaction between a projectile particle i and target nucleus j, as 
described in section [T. 1.31 

The non-resonant form of (o-f ) y , equation i 1.291 contains the quantity S{Q) 
which is the astrophysical S factor at zero energy. Why zero energy? Sim- 
ply put, the temperatures that we find in the centers of stars range from lO'^K 
for hydrogen burning to a few times 10^ A' for helium burning, or roughly 1- 
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20keV. Though this represents the mean thermal energy of the particles, it is 
an order of magnitude or more below the energies available to nuclear physics 
experimentation. In Figure 11.21 we show representative experimental efforts 
to measure the astrophysical S factor for a nonresonant reaction. Clearly, the 
extrapolation to zero energy needs to be done with some care given the exper- 
imental uncertainties alone, added to the analytic approximations. For more 
precision, additional terms allowing for the value of S at stellar energies can 
be derived from the experimental data for some reactions. 

For resonant reactions, equation [L28] will be effective if the resonance lies 
within the range of energies of the interactions (or, more precisely, near the 
Gamow peak) providing that the energy width F, and its position £'i-cs is known. 
In practice, though, we see from Figure [L2] that the relevant energies are poorly 
sampled by experiment in many cases. Thus a "hidden" resonance can affect 
the behavior of S{E) significantly, resulting in poor estimates of {av)ij. Per- 
haps the most important example of the complicating role of resonances in 
nuclear fusion rates, in terms of stellar astrophysics, is the ^^C(a,7)^^0 re- 
action that is mostly responsible for setting the C/O ratio in the Galaxy. There 
are low-energy resonances (at about 10 Me V) that are po orly mapped, and 



affect the cross-section at lower energies; see lEl EidI (120051) for a discussion. 
Depending on the data and extrapolation method used, the cross section at tem- 
peratures relevant for helium burning (a few x IQ^K) can vary by fairly large 
factors. Figure [Ol shows these cross sections for a few formulations of this 
reaction. 

With the above caveats, there are several recent tabulations of nuclear reac- 
tion rates and thei r dependence on temperatu re and density. Still, the "classic" 



reference remains ICaughlan & Fowleij ( 119881) and refer ences therein, with u p 



dated rates available on the NACRE database / website (lAngulo et ami999h . 



"Thermal" neutrino emission 

One of the most important coolants for white dwarfs during their early phases 
is neutrino emission. The processes that produce neutrinos in white dwarfs 
differ from the neutrino production associated with nuclear fusion. In a dense 
plasma, neutrinos can play the role that photons play in more ordinary stellar 
material. When photons are produced in these interactions, they quickly ther- 
malize with the plasma because of electromagnetic scattering. But when neu- 
trinos are produced, they are not thermalized and leave the star, taking their 
share of the interaction energy away with them. Thus neutrino emission is an 
energy sink, rather than a source. 

The ability of neutrinos to act in this way is enhanced when the density 
is high - and so thermal neutrino processes can be important coolants in the 
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Figure 1.3 Variation of th e cross section of t he ^^C(q:,7)^^0 reaction for 
various formulations, from lEl Bid et alj (|2004) . At typical relevant stellar en- 
ergies are approximately a few x 10^ K - indicated by the red line - the range 
of suggested values is significant. See lEl Hid et alJ ( l2004l for details. 



cores of red giants, and in white dwarf stars. The important thermal neutrino 
processes include Bremmstrahlung, where neutrinos are involved in free-free 
scattering instead of photons. "Plasmon neutrinos" are an even more efficient 
energy sink in hot white dwarfs - they result from decay of plasmons within a 
dense plasma. Plasmons are similar to photons, but are coupled with the plasma 
in such a way that they can decay (into a neutrino - antineutrino pair) and still 
conserve momentum and energy. For an excellent a nd compact summar y of 
neutrino production mechanisms in white dwarfs, see lWinget et al.l (120041) . 

The neutrino emission rates are, as is easy to see, not possible to measure 
experimentally. Even if they could be produced under experimental conditions, 
detection of those produced would be an exciting challenge, since thermal pro- 
cesses do not produce coherent neutrino beams. The rates are calculated using 
the Standard Model of particle physics. A current alg orithm for using those 
calculations in the form of energy loss rates in stars is l ltoh et all (Il996!) . Ver- 
ifying those computations does fall into the realm of astrophysics - neutrinos 
are the dominant coolant for most white dwarfs during a significant fraction 
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Figure 1.4 Photon and neutrino luminosity as a function of cooling age for a 
O.615M0 white dwarf with a helium-rich atmosphere. Note that the energy 
loss by neutrinos exceeds the photon luminosity over a significant time in 
this early cooling phase. The plasmon neutrino process dominates Bremm- 
strahlung for this model. 

of their early cooling stages. So, measuring the cooling rate of white dwarfs 
can provide an experimental, if indirect, test of the theoretical thermal neu- 
trino production rates. This measurement can be done in a statistical way by 
looking for features in the luminosity function of a collection of white dwarfs, 
but a more direct way is to measure the cooling rate of a single white dwarf. 
Figure [L2. 21 gives a preview of this, showing that the neutrio luminosity can 
exceed the photon luminosity during a portion of the cooling history of a white 
dwarf. We will discuss this further in a later section, but clearly a sensitive 
asteroseismic probe of white dwarf structure can enable this measurement. 




1.2.3 Opacities: K{p,T,Xi) 

For regions of the star that are not convectively unstable, reliable calculation 
of the opacity can be non-trivial, especially in the outer layers where the ra- 
diative opacity is dominated by atomic processes. The Rosseland mean opac- 
ity (equation 1 1.3 II ) hides these difficulties in the frequency dependent opacity 
K,y. These processes include electron scattering (which is easy as it is not fre- 
quency dependent), free-free scattering, bound-free absorption, and bound- 
bound absorption. At cooler temperatures opacity becomes significant for 
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Figure 1. 5 Frequency-dependent absorption coefficient for iron at logT — 
5.4, from lRogers & IglesiasI ( fl992h 



solar-type stars, and even cooler stars require accurate treatment of molecular 
absorption by CO, OH, H2O, CH4, etc. An issue in computing the Rosse- 
land mean opacity is that the frequency-dependent cross section is an ex- 
tremely rapidly varying function at temperatures near and below complete i on- 
ization. Figure 11.51 illustrates this for iron, from [Rogers & IglesiasI (119921) . at 
2.51 X lO^K, with the relevant atomic transitions as labeled. 

Considering the fact that one needs to include all species that contribute to 
the opacities, integrals for the Rosseland mean opacity must cover dozens of 
elements (each in their appropriate ionization and excitation states). Indeed, 
advances in modeling atomic states and in computational scale resulted in sig- 
nificant revision o f the atomic opacit i es use d in stellar interiors calcu lations in 
the mid 1990s bv lRogers & IglesiasI d 19921) and lSeaton etal.l (Il994l) . It is no- 
table that the need for updated opacities (particularly the contribution of iron) 
was pointed out in part as the result of s eismic modeling of Cepheid variables 
(Simon, 1982 : Iglesias & Rogers , 199l l) as we will see in a later section. 

For cooler stars, the contribution of molecular transitions (and, at lower tem- 
peratures, dust opacity) complicates the calculation further. Figure 11.61 from 
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Figure 1.6 The absorption coefficient for CO as a function of wavelength, 
from Alexander et al. QQQ^. 



Alexander et alj (l2003b . shows a sample run of the absorption coefficient (per 



molecule) for CO. The relative s trengths of the variou s low-temperature opac- 
ity sources are summarized bV iFerguson et alj (l2005b . 



1.2.4 Summary 

Overall, stellar models do a very good job of describing many of the observ- 
ables of stars with the cuiTent state-of-the-art of physics inputs (as long as 
we continue to not mention convection). When asteroseismic efforts are in- 
troduced, though, gaps in our physics understanding are quickly exposed, and 
this in turn can shake our confidence in the way we treat the input physics 
for stellar models. The cuiTent section reviewed the approach taken in modern 
stellar modeling - and in various places, improvements have been made based 
on seismic analysis. 

In the next sections, we will show how observations of pulsating white 
dwarfs can test or constrain the physics of dense matter under conditions when 
crystallization might occur. Those pulsating white dwarfs also place constraints 
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on interesting nuclear cross sections (in particular, for the ^'^C{a, j)^^0 reac- 
tion, and on the neutrino emission rates. Cepheid pulsation systematics de- 
manded recalculation of radiative opacities that include a realistic treatment of 
heavy elements. Observations of pulsating sdB stars probe other effects such 
as radiative levitation (as we will see), and solar-like oscillations of main se- 
quence stars, as well as white dwarfs, provide new probes of convection in stars 
that may finally shake our reliance on the antiquated mixing length theory. 



1.3 Seismology to the rescue: feedback between pulsation 
studies and input physics 

In the previous section we reviewed some of the more difficult or problem- 
atic aspects of the physics of stellar interiors and how that physics can be 
(compactlyO implemented in stellar models. Now we can discuss how aster- 
oseismic studies have impacted this area of stellar astrophysics (mostly for 
the better). The asteroseismic signal of greatest utility is the set of pulsation 
frequencies - using relatively simple adiabatic pulsation theory provides the 
tool we can apply. We'll cover some examples without trying to be exhaustive, 
drawing heavily on results from compact pulsating stars like white dwarfs - 
since that is an area that has been around for a while and one that I'm familiar 



with. lAlthaus et al.l (120101) provide a comprehensive review white dwarf evo- 



lut ion in the conte xt of asteroseismology that is complementary to the review 



bv lKawaled(ll995b 



We start with white dwarf asteroseismology and how we can probe the 
crystallization of stellar interiors, and also how seismic probes of the inter- 
nal composition profile can constrain important nuclear cross sections. Next 
we'll briefly discuss how timing of white dwarf pulsations over many years 
can provide a measurement of the cooling rate of individual stars, and there- 
fore constrain neutrino emission rates. For variety we'll then consider how the 
theory of Cepheid pulsations led to a revision of the radiative opacity calcula- 
tions. 



1.3.1 White dwarf crystallization revealed by asteroseismology 

In Section [1.2. II we told the story of how Salpeter ( 1961 *) concluded that un- 
der conditions found within white dwarf stars, the i nterior could un dergo a 
phase transition and develop a crystalline core. What ISalpeten (1196 ih did not 
anticipate was that we would find that white dwarfs undergo nonradial pulsa- 
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tion when the (pure hydrogen) envelopes reach temperatures where hydrogen 
is partially ionized at depth. 

Generally at Toff ~ H, OOOK to 13, OOOi^ the envelope partial ionization 
drives convection and also pulsations, through processes to be reviewed later 
The historical development of white dwarf asteroseismology, which dates from 
the dis covery of the fi rst pulsating white dwarf by Arlo Landolt in the late 
1960s ( LandolJ, 1968); for an excel lent overview of the current state of the 
subject, see IWinget & Kepler! (120081) and Fontaine & Brassard (l2008l) . There 
is a correlation between the luminosity of a w hite dwarf and the core temper- 
ature that follow s closely the seminal work by iMestell (119521) and updated by 
Van Hornl(ll97lb : 



L 

— = 1.7 X IQ- 

Lq 



M 



(1.47) 



where the core temperature is in units of lO'^K. Recalling Section [l.2.1l and in 
particular equation 1 1.461 the central temperature at crystallization for oxygen 
is approximately 3.4 x lO^K, which is at a luminosity of only sa 2 x IO^^Lq. 
This is at the lower limit of observed white dwarf luminosities, suggesting that 
any white dwarfs that are cool enough to crystalize will be too cool to pulsate. 
However, one must remember that the radius of a white dwarf is determined 
almost entirely by its mass (the famous mass-radius relationship for degenerate 
stars) - so if the mass of a white dwarf at that luminosity is large enough, its 
would still require a sufficiently high T^s to be seen at this luminosity. 

In fact, white dwarfs with masses greater than about IMq can crystalize 
at effective temperatures that place them in the pulsational instability strip. 
This possibility was first expl ored by Mike Montgomery and Don Winget 
(IMont^omerv & Winged 1 19991) following the discovery of a massive p ulsating 



white dwarf, BPM 37093, with a mass of approximately I Mr^ (|Kanaan et al. 

992h . Computations of evolutionary white dwarf models bv lMonteomerv & Winget 
(Il999h (see Figure [LTT i compared with the then-known pulsating white dwarfs 



confirmed the theoretical expectation that BPM 37093 should be a pulsating 
white dwarf with a crystalline core. The seismic diagnostic of a crystalline core 
is the fact that nonradial modes propagate quite differently through solids than 
through gasseous material. For g— modes such as we see in white dwarfs, the 
eigenfunctions carry very small amplitude in crystalline material when com- 
pared to non-crystalline white dw arf cores, which in turn a ffects the pulsation 
frequencies in a measurable way. iMontgomerv & Wingeti (Il999h demonstrate 
how a of g— mode overtones can determine the mean period spacing, and cou- 
pled with other data can reveal whether the core is crystalline or not. 

In an effort to determine the frequencies of as many modes as possible. 
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Figure 1.7 Crystalline core fractions for white dwarf models as a func- 
tion of mass and T cb compared with some pulsating white dwarfs, from 
iMonteomerv & Winget ( 1999) . 



BPM 37093 was the primary target for two Whole Earth Telescope campaigns 
that ultimate ly were able to expose se veral 1=2 q— m odes in the star jKanaan et al 
200dl2005l) .Fi gure ll.81 from Kan aan et alj (l2005h shows the measured period 
spacing superimposed over various models. Modehng of these pulsations con- 
cluded that there is strong evidence that the core of BPM 37093 is indeed crys- 
talized, in accordance with theoretical expectations. With the conclusion that 
BPM 37093 was largely crystalline, the public became highly engaged. This 
asteroseismic result was announced in February 2004, making the discovery 
of an enormous "diamond star" a terrific Valentine's Day story! 



1.3.2 The C(a, 7) °0 reaction rate - constraints from white 

dwarf seismology 

We saw in Section 11.2.21 that the reaction that is principally responsible for 
the C/O ratio in the Universe has large uncertainties associated with it because 
of the difficulty mapping low-energy resonances. The ^^C(a, 7)^^0 reaction 
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Figure 1.8 The observed period spacing in B PM 37093 compared with rele- 
vant stellar models, from lKanaan et al.) ( I2005h . 



occurs at slightly higher temperatures than the triple— a reaction that produces 
^^C during core helium burning, and so the rate of production of ^^O increases 
towards the end of helium core burning. White dwarfs therefore should have 
a decreasing C/0 ratio as one progresses from the outer core to the inner core 
reflecting this evolution. The central C/0 ratio will therefore depend on the 
^^C(a, 7)^^0 reaction rate. 

White dwarf pu l satio n periods principally depend on the envelope structure 
jWinget & Keplerl 12008) . but the core properties also influence the periods. 
In matching observed periods with model periods, near-surface composition 
gradients (accentuated by diffusion processes) produce mode trapping, allows 
determination of the surface layering structure. But composition gradients in 
the core (c aused by the increasing production of ^"^O, for example) are also 
influential jMetcalfe et all boosi; [Montgomery et aLiliooj) . 

With several DB (helium atmosphere) white dwarfs having successful ground 
based observing campaigns that revealed a large number of modes, several 
groups have investigated whether seismic probes could in fact measure the core 
oxygen abundance and the composition gradients in the core caused by the se- 
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Figu re 1.9 Oxygen abundance profile in a pulsating white dwarf interior, 
from lMetcalfd OOOSi) . The wide (gray) shaded regi on is the range of pro files 
given the ±la range from the NACRE compilation dAngulo et allll99^ : the 
seismic constraints are given in the narrower in band for two white dwarf 
pulsators. 



quential stages of helium burning. iMetcalfe et alJ (120011 120021) showed that by 
varying the assumed ^^C(a,7)^^0 reaction rate, the core oxygen mass frac- 
tion in th e remnant wh i te dw arfs could chang e significantly. Using data from 



GD 35 8 (IWinget et al.L 1 1994 . and CBS 1 14 (iHandleret al.L l2002h . Metcalfe 



(120031) concluded that the rate for this proble matic reaction are in fact in line 
with recent determinations discussed by .El Eid et al.l (120041) . The seismic re- 
sults are shown in Figure [L9l 



1.3.3 Measuring white dwarf cooling via asteroseismology 

In section 11.2.21 and in particular Figure 11.2.21 we saw that the cooling (and 
therefore overall evolution) of white dwarfs can be dominated by plasmon 
neutrino emission during the drop in luminosity from the planetary nebula nu- 
cleus phase into the low-luminosity phase. By measuring the cooling of an 
individual white dwarf, whose global properties are otherwise constrained by 
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spectroscopy and asteroseismology, we can in principle test the theoretically 
expected rate via observation. 

Measuring cooling rates via pulsations 

Consider the accuracy to which we can measure most stellar properties through 
spectroscopy. On a good day, we can measure, for example. Toff to 0. 1 %. With 
knowledge of the distance to a star, we might gain a precision in luminosity of 
perhaps 1%. Now, if we want to measure changes in those quantities brought 
about by the slow secular evolution of the star, which occurs on time scales of 
milhons or bilhons of years, we'd have to stick around for a thousand years or 
more before being able to detect any changes at all. 

But through time-series photometry or spectroscopy, we can measure os- 
cillation frequencies to a much higher precision, especially in those stars that 
show self-excited (and therefore extremely coherent) pulsations. Stars for which 
we can get parts-per-million accuracy in periods include the pulsating white 
dwarfs, 5 Scuti stars, and other classical pulsators such as RR Lyr and Cepheid 
pulsators. With measurement of a global property of a star (a pulsation period) 
at that level of accuracy, we can indeed hope to measure a change in the pul- 
sation period, caused by secular evolution, in a relatively short period of time 
(within the career of a graduate student, for example). 

The direct measurement of a pulsation period change requires patience and 
a cooperative star that shows tractable phase or ampUtude variations over many 
years. While one could in principle measure the pulsation period itself directly 
at two widely-spaced epochs, a more integrated method would be to look for 
continuous phase changes caused by the increasing or decreasing period. If we 
measure the period very precisely, and use it to predict the phase in the future 
under the assumption that the star's period is not changing with time, then any 
real changes in the period within the star will cause the phase to advance (or 
retreat) in a way that will reveal the period change. The traditional method for 
accessing these phase changes is through a so-called {O — C) analysis, where 
(O — C) is the difference between the observed phase at a given time, and that 
calculated assuming that the pulsation period is constant. 

The phase, as determined from the time of next maximum in the Ught curve 
at a given epoch, is given by 

dP 

imax.i = imax.O + « X Po + i^J'-J- (1-48) 

(It 

for the ith maximum after an agreed-upon first maximum, where we have 
expanded the equation through a Taylor expansion assuming dP/ dt is small. 
If the period itself does not change with time, the predicted maximum of the 
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Figure 1.10 The change in (O — C) as a function of time for constant rate of 
increase in period (upper curve) or decrease in period (lower curve). 



ith cycle will be perfectly accurate with only the first two terms in the above 
equation, and the value of (O — C)i will be a constant at zero. If there is an 
error in the period, then (O— C) will increase linearly with time (if the assumed 
period is too short) or decrease linearly with time (if too long). 

If the Pq is not in error but the period is changing (slowly) with time, then 
equation !! .481 vields 

1 rlP 

iO-C) = {Atf-— , (1.49) 

where At is the time elapsed between the chosen epoch and the measurement 
of the current time of maximum imax- Clearly, the departure of (O — C) from 
zero increases as (time)^ goes by, enabling quadratically increasing sensitivity 
to dP/ dt over time. A change in (O — C) will have a parabolic form if dP/ dt 
is constant, as shown in Figure [TTTOl 

Since the pulsation period is some measure of the global structure of the star 
(such as the radius as it affects the mean density), the period of pulsating stars 
much change as they evolve, and must change on approximately the evolution- 
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ary time scale. Assuming a a time scale for period change (inversely propor- 
tional to P / P) of 10^ years, and a precision of measuring <,„ax of 10 seconds 
(both typical for young, hot white dwarfs), one finds that the time needed to 
detect an evolutionary change is only about two years for a Ict result, and 3 
years for a 3a result. 

Of course, these expectations are for ideal circumstances. For the above pro- 
cess to work, we need to ensure that the frequencies are precisely known, and 
that the analysis is not hampered by cycle-counting errors during the inevitable 
daily, monthly, and seasonal gaps in the data. In addition, not all oscillation 
modes are not as stable in phase and/or amplitude as needed to reveal this 
effect. Nonlinear interactions between modes, shorter-timescale phenomena 
such as convection or magnetic fields, and other effects can hamper our ability 
to do the necessary high-precision cycle-counting required. All of these ef- 
fects would only lead to larger values for P, so that the estimates for the rate 
of period change (and therefore the evolutionary time scale) are always lower 
limits to the true time scales of evolution. 

White dwarfs evolve at roughly constant radius, so their evolution is a cool- 
ing process, with the fading over time resulting entirely from a decreasing 
core (and envelope) temperature. The pulsations are modes, which are sen- 
sitive to the Briint-Vaisala frequency. That, in turn, is set by the temperature 
stratification in the star - hence the cooling of the star leads to a decrease in 
the Briint-Vaisala frequency, and an increase in the pulsation period. Thus for 
white dwarfs we expect to see periods increase with time, on time scales equal 
to the cooling rates of those stars. 

An excellent demonstration of this is the pulsating DA white dwarf G117- 
B15A. This star has been continually monitored for well over 30 years, mostly 
by the Brazilian astronomer S. O. Kepler (no relation). Figure [TTTT] shows the 
(O — C) diagram for this star, and the (now) statistically s ignificant measure- 



ment of the rate of cooling of this star dKepler et aUKOOSh . This star is signif- 



icantly cooler than the example cited above, and has an expected cooling time 
scale of approximately 10^ years. 

Application to wliite dwarf pulsators 

There are two classes of pulsating white dwarf that overlap with this neutrino 
cooling phase. At the higher-luminosity end are the GW Vir stars, with the 
prototype PG 1159-035. The luminosities of those stars are generally higher 
than of interest for neutrino cooling studies, but th e lower-luminosity m em- 



bers are in the range where L^, > L^. In particular. lO'Brien et alJ (Il998h and 



O'Brien & Kawaled (l2000l) show that the GW Vir star PG 0122H-200 would 



be a good candidate for measuring the cooUng rate and therefore constrain- 




Figure 1.11 An (O — C) diagram for the 215s g— mode in the pulsating DA 
white dwarf G117-B15A, revealing the upward parabol a that results from a 
positive rate of period change, from lKepler et alJ ( l2005h . The time scale for 
period change is approximately 1.7x10 years. 
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Figure 1.12 Relative photon and neutri no luminositie s in th e region where 
we find DB pulsating white dwarfs, from lWinget et al.l ( 120041) . 



ing neutrino emission. Unfortunately, a recent {O — C) analysis suggests that 
the period of PG 012 2+200 is changing much faster than the models predict 
(IVauclair et al.L l201 ih . suggesting that other factors are influencing the pulsa- 
tions. 



Winget et al.l (120041) show an example of this approach (see Figure 11.131 1. 



and demonstrate that the rate of period change can be an effective probe, since 




Figure 1.13 Rate of period change in various DB white dwarf models, from 
IWinget et alj ( 120041) . The solid curve connects modes in a model with the 
standard neutrino emission, while the dotted line is for modes in a model 
with neutrino emission suppressed. 



the rate of period change depends nearly linearly on the neutrino rates used in 
the calculations. 

More promising are stars at the lower end of the range of luminosities where 
neutrino emission dominates photon cooling (see Figurel 1.2.211. E xpanding the 
lower luminosity side in Figure 11.121 (from Winget et al. ( 2004!)) we see that 
two well-studied pulsating DB white dwarf stars straddle the luminosity where 
neutrino emission yields to photon emission as the dominant cooling process. 

In particular, EC 20058 should be experiencing neutrino cooling that ex- 
ceeds photon cooling by nearly a factor of two. The influence of neutrinos on 
the cooling rate can be estimated by computing evolutionary models of the 
stars with the current neutrino rates, and with modified neutrino rates. Self- 
consistent calculations with and without neutrino emission, for example, both 
constrained to match the observed pulsation periods themselves, can then be 
probed to provide the expected rate of period change as a function of the neu- 
trino emission rate. 



Sulhvan et al.l (l2008h report on Whole Earth Telescope observations of EC 



20058, which is part of a long-term observing campaign that spans the time 
from its discovery as a pulsator in the mid 1990s through the present. The star 
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is an extremely stable pulsator (in terms of phase and amplitude changes) and 
should, in the near future, provide the first strong constraints on the rate of 
plasmon neutrino emission from a dense plasma. As Denis Sullivan would say, 
"stay tuned!" 



1.3.4 Cepheid masses as opacity probes 

Here, I will briefly recount a saga from the annals of stellar pulsation the- 
ory that was an early demonstration of the power of asteroseismology (though 
it was yet to be named that) to expose and fix problems with stellar interior 
physics. The centerpieces of this story are multimode pulsators: the so-called 
"beat Cephieds" and "bump Cepheids." The beat Cepheids (also sometimes 
called double-mode Cepheids) pulsate in the (radial) fundamental and first- 
overtone modes. The period range of the beat Cepheids is in the 2-7 day range, 
and the ratio of the period of the first overtone to the fundamental (Pio) ranges 
from about 0.695 (at the long period end) to 0.715 at the short period end. 
The "bump Cepheids" display a secondary maximum that progresses slowly 
in phase with respect to the fundamental pulsation mode; the bump is the man- 
ifestation of a near-resonance between the pulsation of the second overtone P2 
and the fundamental Pq. The bump Cepheids are a more homogeneous class 
with periods close to 10 days, and P20 very nearly 0.5 - that is, close to a 1:2 
resonance. 

As summarized very nicely in iMoskalik et al.l ( Il992h and ISimonI ( Il987h . 
evolutionary models of beat Cepheids with the proper fundamental period Pq 
had period ratios Piq that were significantly larger than the observed value. 
To reach the observed Piq models using the then-standard opacities (i.e. in the 
early 1980s) required unrealistically low masses; even so, the models had a 
steeper dependence of Piq on Pq than what is observed. For bump Cepheids, 
the period of 10 days corresponded to a mass of 6.3Af0 - but then the theoret- 
ical value of P20 was too high to explain the observed resonance. Thus as of 
1981, the evolutionary model masses for Cepheids did not match the pulsation 
masses, with the discrepancy being large enough (a factors of two more) to 
make this a "famous problem." 

In 198 1 and 1982, Norman Simon made a bold statement that these problems 
with Cepheid masses could be solved if the opacity caused by heavy elements 
was being underesti mated by a fac t or of 2 -3 in the then-current generation of 



opacity calculations (iSimonL 119811 Il982h . He found that by making such an 



increase, the resulting models were able to explain the observed periods and 
period ratios with the same mass as required by evolutionary calculations. This 
same augmentation also would provide a driving mechanism for the pulsating 
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Figure 1.14 Interior opacit y in a Cepheid model from the Los Alamos compi- 
lation ( ICox&Tabo3.ll976h . along with the suggested enhance ment needed t o 
reconcile evolutionary and pulsation models of Cepheids, from lSimoi] ( ll982h . 




B stars known as /3 Cepheids, another "famous problem" of the era. His 1982 
paper, "A Plea for Reexamining Heavy Element Opacities in Stars" states 

It is thus quite possible that, by the single stroke of augmenting the heavy element 
opacities by factors of 2-3, we can bring into line with the theory of stellar structure 
and evolution not only the double-mode and bump Cepheids, but the /3 Cepheid 
pulsators as well. 



Simon 's augmented opacities from 1982 are shown in Figure [l~T4l from lSimon 
(Il982h . The enhancements at temperatures from lO'' * and above, he postu- 



lated, should result from a more careful treatment of b-f and b-b transitions in 
heavy elements than were implemented in the standard opacity tables available 
at that time. The augmented opacities in Figure [TTl4]pi"oduced models that fit 
the beat and bump Cepheid phenomena, and solved other problems in stellar 
pulsation theory. 

This suggestion, in part, motivated Carlos Iglesias and Forrest Rogers, of 
Lawrence Livermore National Laboratories, to undertake a new computation of 
astrophysical opacities. They eventually turned the computational machinery 
of their laboratory onto the problem, and found that indeed, just as predicted 
by Simon nearly a decade earlier that inclusion of an improved treatment of the 
physics of bound st ates in heavy e lements led to an enhanced opacity. Figure 
ll.lSl shows that, as Simon ( 1982 1) predicted, the opacity enhancement was a 



factor of 2-3 at the temperatures in that e arlier work. 



Using the (then) new OPAL opacities, iMoskalik et al.l (Il992h revisited pul 
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Figure 1.15 Earl y OPAL opacities relevant to Cepheid interiors, compared 
to analytic fits bv | Stelling werf ( 1974, 19 7^ to the L os Alamos compilation 
iCox & Taboiill976l) . from ilglesias & Rogeri jl991h . 



sation and evolution models of Cepheids (and (3 Cepheids) to verify that the 
OPAL opacities did indeed solve the Cepheid mass problem. Their results 
provided a dramatic closure to the problem, confirming that the newer in- 
put physics resulted in a much closer match between the masses of Cepheids 
from evolutionary consideration with the masses determined through pulsa- 
tion period and period ratio matching. Figure [TTSl for example, shows tracks 
for Cepheid models with the Los Alamos opacities (dashed lines) falling well 
above the observed beat Cepheid period ratios in the P10--P0 plane, while mod- 
els with the OPAL opa cities fit the data quite closely and at reasonable masses. 



Moskalik et alJ (Il992h also show concordance for the bump Cepheids. 
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Figure 1.16 Period ratios for bea t Cepheid models compared with observa- 
tions, from lMoskalik et al] ( Il992h . Dashed lines correspond to models with 
older opacities, solid lines are tracks for models with OPAL opacities. 



1.4 Summary 



The previous section gave only a few illustration of the ways that stellar pulsa- 
tion and asteroseismology have been able to teach us some basic physics. Many 
other examples are in the current literature and form the core of active research 
in asteroseismology. We haven't even mentioned the impact that helioseis- 
mology has had on the d i scovery and characte rization of neutrino oscillations 
jBahcall & Ukichl Il988t iBahcaU et all I2OO2I) o r on probes of the details of 



convec tion in the outer layers of the Sun such as in lChristensen-Dalsgaard et al 



(1199 ll) . Details of convective efficiency in w hite dwarfs, too, can be exposed 

throu gh analysis of the shape of the pulsations jMontgomervil2007tlMontgomerv et al, 
2Q10l) . Even more fun is the possibility that measuring the evolution rate of 
white dwarfs could reveal the pr esence of axions or other exotic particles, as 
suggested bv lBischoff-Kim et al.i (.2008.) and others. 

These results are meaningful for our selfish purposes of improving stellar 
models for application to other areas of astrophysics, but also have value well 
outside the core astrophysics disciplines - atomic physics, nuclear physics, and 
condensed matter physics in particular. All this from a handful of coefficients 
in four differential equations ! 
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